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In this talk, we examine the existing calculations of QCD evolution kernels for the 
scale dependence of two sets of twist-3 quark-gluon correlation functions, T q p(x,x) and 
T (<T j(x,x), which are the first transverse-momentum-moment of the naivc-timc-revcrsal- 
C ~ ) odd Sivers and Boer-Mulders function, respectively. The evolution kernels at the leading 

r- H order in strong coupling constant a a were derived by several groups with apparent dif- 

ferences. We identify the sources of discrepancies and are able to reconcile the results 
from various groups. 

(N 

> 
X 

Quantum Chromodynamics (QCD), as a fundamental theory of strong interaction 
dynamics, has been very successful in both lattice QCD calculations of static hadron 
properties and perturbative calculations of short-distance dynamics observed in 
high energy scattering processes. In particular, with the aid of QCD factorization 
theorerrP, QCD perturbation theory has done an excellent job in interpreting data 
from high energy scattering with identified hadrons, whose nonperturbative dy- 
namics at the hadronic scale (~ 1 /fm) are represented by process independent and 
well-determined parton distribution and fragmentation functions^ - ^. The predictive 
power of QCD perturbation theory for cross sections with identified hadrons is an 
immediate consequence of QCD factorization and our abilities to predict the vari- 
ation of these nonperturbative but universal distribution and fragmentation func- 
tions when they are probed at different hard momentum scales. The success of QCD 
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perturbation theory provides the well-controlled and calibrated sub-femtometer or 
even attometer "scope" to probe partonic structure and dynamics inside a hadron 
of bounded quarks and gluons. 

On the other hand, the large size of observed transverse single-spin asymmetries 
(SSAs), An = (ct(st) — cr(— st))/(c(st)+o'(~st)), defined as the ratio of the differ- 
ence and the sum of the cross sections when the spin of one of the identified hadron 
st is flipped, came as a surprise, and had posed a challenge for researchers in this 
field for some time^. From the parity and time-reversal invariance of the strong in- 
teraction dynamics, the measured large asymmetries in high energy collisions should 
be directly connected to the transverse motion of partons inside a polarized hadron. 
Transverse spin physics has attracted tremendous attention from both experimen- 
tal and theoretical communities in recent years^. Experimental measurements of 
the asymmetries and the investigation to understand the underlying dynamics have 
provided and will continue to provide us new opportunities to explore QCD and the 
hadron structure far beyond what we have been able to achieve. 

Two complementary QCD-based approaches have been proposed to analyze the 
physics behind the measured asymmetries: the transverse momentum dependent 
(TMD) factorization approactP~OS and the collinear factorization approach 110 ! 11 ! 12 ! 
In the TMD factorization approach, the asymmetry was attributed to the spin 
and transverse momentum correlation between the identified hadron and the ac- 
tive parton, which are represented by the TMD parton distribution or fragmenta- 
tion function. On the other hand, in the collinear factorization approach, all ac- 
tive partons' transverse momenta are integrated into the collinear distributions, 
and the explicit spin-transverse momentum correlation in the TMD approach is 
now included into the high twist collinear parton distributions or fragmentation 
functions. The asymmetry in the collinear factorization approach is represented by 
twist-3 collinear parton distributions or fragmentation functions, which have no 
probability interpretation, and could be interpreted as the quantum interference 
between a collinear active quark (or gluon) state in the scattering amplitude and a 
collinear quark (gluon)-gluon composite state in its complex conjugate amplitude. 
The relevant TMDs and the twist-3 quark-gluon correlation functions, while they 
are both nonperturbative, are closely related to each other. In general, the collinear 
twist-3 correlation functions are proportional to the parton transverse momentum 
/cj^-moments of TMDs up to the uncertainty of ultraviolet renormalization of com- 
posite operators defining the moments of TMDs. The TMD factorization approach 
is more suitable for evaluating the SSAs of scattering processes with two very dif- 
ferent momentum transfers, Qi ^> Q2 > Aqcd, where the Q2 is sensitive to the 
active parton's transverse momentum, while the collinear factorization approach is 
more relevant to the SSAs of scattering cross sections with all observed momentum 
transfers hard and comparable: Qi ~ Q 3> Aqcd- Although the two approaches 
each have their own kinematic domain of validity, they are consistent with each 
other in the regime where they both appljUHEH 

Both factorization approaches necessarily introduce a factorization scale, fi 3> 
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Aqcd, to separate the calculable short-distance perturbative dynamics from the 
long-distance nonperturbative physics of the observed cross sections or the asym- 
metries. Since the physical observables, the cross sections or the asymmetries, are 
independent of the choice of the factorization scale, the scale dependence of the 
nonperturbative distributions^'"'^', either TMD distributions or twist-3 collinear 
distributions, must match the scale dependence of corresponding perturbative hard 
parts. That is, the factorization scale dependence of the nonperturbative distribu- 
tions is perturbatively calculable and is a prediction of QCD perturbation theory 
when /j, 3> Aqcd- For example, the scale dependence of the leading power parton 
distributions obeys DGLAP evolution equations whose evolution kernels are per- 
turbatively calculable, and has been very successfully tested when the scale varies 
from a few GeV to the hundreds of GeV. 

In this talk, we present a general evolution equation of the twist-3 quark-gluon 
correlation functions that are responsible for the SSAs in the collinear factorization 
approach^, and a detailed discusion on the evolution of two twist-3 quark-gluon 
correlation functions T q p(x,x) and T^" F (x, x). These two correlation functions are 
defined as EMI 

T q Ax,x) = J d ^^e" p+ y^(P lST \^ q (0) 1 + [e s - ann F+(y^)] ^)\P,s T )(,l) 
T$( X) x)= [ ^^ e - p+ y^-J2^ s T\M0) [a<*+F a +(y;)]Tp q (y;)\P,s T )(,2) 

where the gauge links between field operators are suppressed and e 0123 = 1 is 
used. These two correlation functions are equal to the first fc^-moment of the two 
well-known naive-time-reversal-odd TMDs, the Sivers function /^(x, k^J^ and 
the Boer-Mulders function hj^(x, kjj*\ respectively. The scale dependence of these 
two twist-3 quark-gluon correlation functions has been studied recently by several 
groupfOQ 21 l 22 l 23 l However, there are discrepancies between these results, particularly 
for the evolution of T q ^p(x,x) (also often refer to as Efremov-Teryaev-Qiu-Sterman 
(ETQS) function). We show in this tallP^that these discrepancies can be resolved, 
and in addition we also present the evolution equations for the quark-gluon corre- 
lation function T q p(x, x). 

2. Collinear factorization approach to SSAs 

In this section, we give a general discussion of collinear factorization approach to 
perturbative QCD treatment of SSAs of cross sections with one large momentum 
transfer Q. 

A scattering cross section on a hadron is proportional to the square of the 
hadron's scattering amplitude, which sums over all partonic scattering amplitudes 
that share the same initial- and final-states but with different number of active 
partons, as illustrated in Fig. [T] With one large momentum transfer Q, the hard 
scattering is localized to a distance scale of 1/Q. Since pulling an extra physically 
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Fig. 1. A schematic illustration of the source of contribution to the scattering cross section 
between a simple particle and a hadron of momentum p and spin s. 

polarized partem into the localized collision point is suppressed by the power of 1/Q, 
the cross section for a hadron A to scatter off a hadron B with a large momentum 
transfer Q can be expanded in a power series in 1/Q, 



where Q 2 ~ (k 2 ), (k 2 -,) represents a characteristic scale of the power corrections. 
According to the QCD collinear factorization theorem^, the leading power contri- 
bution to the cross section of the hadron-hadron collisions in Eq. ([3]) is given by 
the square of the scattering amplitude with one active collinear parton from each 
colliding hadron (plus any number of collinear and longitudinally polarized gluons 
responsible for the gauge links), like the first term of the scattering amplitude in 
Fig. [1} The contribut ion can be factorized into a convolution of a localized and per- 
turbatively calculable hard part from the collision between partons a and b, 
and the universal twist-2 collinear parton distribution functions (PDFs), f a /A (and 
to find a parton of flavor a (and b) from the hadron A (and B) 7 as indicated in 
Eq. Q. The leading power contribution in Eq. ([3| contributes to the cross section, 
but, not to the SSA. This is because the SSA is a naively time-reversal odd observ- 
able, we need, in order to generate the SSA, a phase and a spin-flip at the partonic 
scattering. At the leading power in the collinear factorization, the phase could only 
be generated by the interference between the tree and the loop diagram, and the 
spin-flip could be achieved by the quark mass. That is, the asymmetry generated 
in this way must be proportional to a s rriQ /px, which is a power suppressed small 
number, and is not sufficient to explain the data on the SSA^. 

In the QCD collinear factorization approach to SSAs, the asymmetry is gen- 
erated by the quantum interference between the first and the second scattering 
amplitude in Fig. [TJ which is the first power correction to the spin-dependent cross 
section and is capable of generating the right phase and required spin flip between 
the single parton and the two-parton composite state to generate the SSAl 10 ! 11 ! It 
was argued in Ref. [M] that such a power correction to the single jet/particle in- 
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elusive cross section of A(pa,s±) + B(jps) — > h(p) + X can be factorized in the 
same way as the leading power term, except that the PDF of the polarized hadron 
is replaced by a twist-3 quark-gluon correlation function Tupy^s) with a = q,q, 
or g as in Eq. (|5j), or a twist- 2 transversity distribution h a /A{s) with a — q, or q to 
take care of the spin-flip while the other PDF is replaced by a spin-averaged twist-3 
quark-gluon correlation function T^ F y B to take care of the phase generation (or 
the twist-2 fragmentation function is replaced by corresponding twist-3 collinear 
fragmentation function 2 -^ in the case of inclusive single hadron production, which 
will not be discussed in this tahf^). 

The twist-3 distributions T a< p(x, x) and T^p(x, x) represent the long-distance ef- 
fect of the quantum interference between a scattering amplitude with a single active 
parton and the one with an active two-parton composite state, as shown in Fig. [I] 
Although the quantum interference between the partonic scattering amplitudes with 
different number of active partons is suppressed by a power of 1/Q, its contribution 
to the SSA could be enhanced in certain parts of the phase space, in particular, 
in the forward region of the polarized hadron, which is a natural feature of twist-3 
contributions^. The predictive power of the approach relies on the universality and 
our knowledge of the twist-3 quark-gluon correlation functions, which are defined 
in terms of matrix elements of twist-3 operators' 11 ! 18 !, e.g., those in Eq. pi). Like the 
usual twist-2 PDFs, these twist-3 parton distributions relevant to the transverse 
SSAs are non-perturbative and universal, and should be extracted from experimen- 
tal data on the SSAs. However, unlike the usual twist-2 PDFs, these distributions 
do not have the probability interpretation and are not necessarily positive. On the 
other hand, we could interpret these distributions (or the matrix elements) as the 
expectation values of the inserted field operators in Eq. ([2]), which corresponds to 
the color Lorentz force and magnetic force experienced by the active partorP3. 



3. Evolution and evolution kernels 

The evolution or the factorization scale dependence of these twist-3 parton corre- 
lation functions is an immediate consequence of the QCD factorization formalism 
for physical observables. If one writes the spin-dependent cross section from the 
factorized formula in Eq. ^ or (J6I), schematically, as 

Aa(Q, s ± ) = (l/Q)H NLP (Q/v, a.) ® /(/i) <8> 7» + 0(1/Q 2 ) , (7) 

with the factorization scale /x, one can derive from d/dln(/i)[Aer(Q, sj_)] = the 
leading order evolution equation for a generic twist-3 correlation function T, 

where J ff NLP ( 1 ) and P 2 (1) are the leading order partonic hard part and the evolution 
kernel of the twist-2 PDFs, respectively. Both terms in the brackets on the right 
of Eq. (pi are perturbatively calculable and give the leading order evolution kernel. 
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Since the distributions are universal, the perturbative evolution kernel could be 
derived in many different ways, although the results should be the samdHMEH 



'V' 




Fig. 2. A diagrammatic illustration for the evolution equation of the twist-3 quark-gluon corre- 
lation function. 



In Ref . [TB] , we presented a derivation of the evolution equations of twist-3 quark- 
gluon correlation functions relevant to SSAs from the perturbative variation of these 
correlation functions. The evolution equations can be schematically represented by 
the plot in Fig. [2] The evolution kernels can be perturbatively calculated from 
the first Feynman diagram on the right-hand-side of the diagrammatic equation 
in Fig. [2] with the bottom parton lines contracted by the cut vertices and the top 
parton lines contracted with the proper projection operators derived in Ref. jTHj. For 
example, for the ETQS function T qtF (xi, X2), the cut vertex V 9i f and corresponding 
projection operator V q _ F in the light-cone (LC) gauge are given by 
, + 



V, 



(LC) 



V 



,(LC) 



2P+ 



-7-P 
2 ' 



hi 

pi 



X 2 5 \X2 



P+ 



(ie*^™) [-g all ] C q . 



(ie STpnTl ) C, 



<i > 



(9) 
(10) 



where (C,)° = (t% and (C,)£ = 2/(N 2 - l)(f 
of quark and gluon fields, and the £ = p + jP + and £2 
fractions of the parent partons, as shown in Fig. [2] 



define the color contractions 
P2 I ' P + are momentum 






fP+P, 




fP+P 2 





Fig. 3. Feynman diagrams contribute to the LO evolution kernel of quark-gluon correlation func- 
tions. 



By calculating the leading order Feynman diagrams in Fig. [3] with the cut vertex 
and projection operator in Eqs. (9) and (10), respectively, we obtained finite contri- 
butions to the leading order evolution kernel and the following evolution equation^, 



dT 9tF (x,x,n) a s f 1 d£ 

= 2tt J 9q q ' F 



(11) 



2 



(1 + z)T qtF (£, x, m) ~ (1 + z 2 )T q Ati, g, M) 
1 - z 



T 



Aq,F 
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where z = x/£ and P qq (z) is the splitting kernel for unpolarized quark distribution 
function, and the quark-gluon correlation function Ta ?i f(^i, %2i fj) is given by^l 



dVi dy~ 2 
Air 



Jx 1 P+y 1 e i(x 2 -x 1 )P+y 2 



(12) 



x(P, ST |V^(0) 7 +7 5 [i8%.F+fo)] ip q ( V 7-)\P,s T ) 



The results derived in Refs. [211 E2] are consistent with ours. But, the evolution 
equation derived later by Braun, Manashov, and Pirnay in Ref. [23 is slightly 
different, 



dT q j?(x, x, /j) 
din fi 2 



2?r 

Nc 
' 2 



y ^Pqq(z)T q , F (£,,£,V.) 

(1 + z)T qiF (£, x, n) - (1 + z 2 )T q , F (£, & /i) 



1-2 

-iV c <S(l - z)T qtF (x,x,n) 



TAq, F (x,£,[J,) 



(13) 



Comparing Eqs. (11) and (13), it is clear that two results differ by two contri- 



butions listed in the third and fourth line in Eq. (13). In addition, there is a sign 
difference in front of the ?a 9 ,f distribution in the second line. This sign difference 
is due to a fact that two groups used a different sign convention for anti-symmetric 
tensor ^ va f i ; we chose e 0123 = 1, while Braun-Manashov-Pirnay used £0123 = 1 
implying e 0123 = — 1. We also noticed that Ref. [27] used the same convention as 
that in our paper, thus they obtain the same sign for the Ta<z,f term. 




SP+P, 



k+k 7 



Fig. 4. Feynman diagrams contribute to the evolution from the interference of a gluon and a 
quark-antiquark state (left), and the usual interference of a quark and a quark-gluon state (right). 



In the fourth line in Eq. (13), the term oc 1/2N C comes from the two Feynman 
diagrams in Fig. |4j The left diagram in Fig. [4] corresponds to the interference be- 
tween a gluon and a quark-antiquark pair. This diagram was missed in our original 
calculation 18 . After including this diagram, we obtain one half of the l/2N c term 
in Eq. (pi 



dT q , F (x,x,n) 



din /1 2 



2tt 



1 d ^ 1 
J2N C 



[(1 - 2z)T q , F {z 



+T Ag>F (z,z-£,/i)], (14) 
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where again the Ta 9 ,f term has an overall sign difference due to our convention 
for e^"/ 3 . The right diagram in Fig. [4] is actually Fig. 7(h) in our original paper 
Ref. [H]. This diagram vanishes if the quark on the left of the cut has a positive 
momentum p + = £P + > 0, which was assumed in the original papei!^. However, 
the £ does not have to be larger than as long as £ + £2 > 0. By calculating the 
contribution from the region where £ < 0, we find that it gives exactly the other 



half of the 1/2N C term in Eq. (13), and reproduce the fourth line of Eq. (13) by 



adding contributions from both diagrams in Fig. [4j 

The term in the third line in Eq. ( 13 ), —N c T q ^ F (x, x), was missed in our original 
paper^. The error was caused by a subtlety in taking the limit x-i — > when we 
evaluate the integration J dx25{x2)x2F(x2) = lirn r2 _>.o £2-^(22) to get the gluonic- 
pole matrix element. The limit, lim^^o X 2F (2^2), would vanish if the function F{x2) 
is finite as X2 — > 0, which is unfortunately not always true in our calculation. We 
find that Fig. |3jb) and (c) have additional contributions to the evolution as^, 



2 

Nc 
2 



d&T q , F (x,x + £2,tJ-) 



X2 
X 2 



X2 



X+X2 — I 



d£,2T q , F {x + x 2 -£ 2 ,x + x 2 ,(J,) 



and 



X2 

e 2 



(15) 
(16) 



respectively. Here the "• • • " includes any regular terms which vanish safely when 
we take x 2 — > 0. The subtlety is caused by the fact that the integration over d£ 2 in 



Eqs. (15) and (16) is singular as X2 — > 0. To evaluate the d^2 integration, we first 
expand T qF (x, x + £ 2 > A*) around £2 = x 2 hi Eq. ( 15 ), 



d 

T qvF (x,x + £ 2 ) = T qiF (x,x + X 2 ) + —Tq tF (%,X + £2)k 2 ^xA& ~X2)-\ -(17) 

C?2 



The integration in Eq. (15) with the first term of the expansion in Eq. (17) gives 



d&T q , F (x,x + X2,Li) 



N c 

= -Y T q,F{X 1 X + X2,tJ-) 



X2 

l-x 



N 



Tq t p(x,X + X 2 , fj) 



X2 
6 



- 1 



(18) 



which goes to — ifT q ^ F (x, x, /1) at the limit X2 — > 0. If one assumes T q ^ F (x, x + £ 2 , A 4 ) 
is a smooth (regular) function, we find that the higher order terms in the expansion 
in Eq. ( 17 ) do not contribute to the evolution in Eq. ( 15 1 when X2 — > 0. So Fig. (3^b) 



gives us an additional contribution — -^Tq iF (x,x,fj,). Similarly we find exactly the 



same contribution from Eq. (16). Adding them together, we have 



dTg iF (x,x,iJ,) 



Fig. [3];b+c)-additional 



-N c T qF {x,x,n), 



(19) 



which reproduces the third line in Eq. ( 13 ) and is what was missed in our original 
paper-^. 
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We now have a complete agreement with the Braun-Manashov-Pirnay result. 
Similarly, our results for flavor singlet evolution are also now consistent with the 
Braun-Manashov-Pirnay result. 



4. Evolution kernel of T^l 

Using the same techniques, we could also derive the evolution equation for the 
other twist-three quark-gluon correlation function T^*p(x, x, //). The calculation is 
straightforward, and the result is 



2tt 



£|{A T P M (z)T$(£,£, M ) 



No 
' 2 



2Tl$(t,x,ri-2zT£J.(£,£,ri 



l q,F 



(<0/ 



q,F\ 



(20) 



N c 6(1 - z) T$(x, x^) + — [2(1 - z)T$(x, x - /x) 



2N r _ 



where AtP m (z) is the splitting kernel for the quark transversity given by 



^TP qq {z) 



2z 



(1-z). 



(21) 



This evolution e qua tion was first derived in Ref. |21j . which contains only the first 



two lines in Eq. (20). The first term in the third line, — N c T^°p(x, x, /x), has exactly 



the same origin as those in Eqs. ( 15 1 and ( 16 1 from calculating diagrams in Figs. 3 



q,F\ 



and (c) with a caution of taking the limit X2 — > 0. The second term, oc 1/2N C , is 
again due to the fact that the Feynman diagrams in Fig. [4] were not included in the 
calculation of Ref. |2T] . 



5. Summary 

We have rederived the evolution equations for both T Qj p(x, x, fi) and T^" F (x, x, t/ipQ 
We resolved the discrepancies in the literature for the evolution of ETQS function 
T q! p(x, x, /i). We understand that such discrepancies were also resolved recently by 
the other two group£32II through careful reexaminations of their original deriva- 
tions in Refs. [3TJ[55], also in Ref. [30] from a different approach. Using the same 
techniques developed in the current paper, we updated the calculation for the evo- 
lution of T^p(x, x, /i) and found two additional contributions which are missing in 
the literature. These results will have important consequences, e.g., in the study of 
QCD resummation for the spin-dependent observable*^. 
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